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O -1 Abstract 

U 

We consider tubular neighborhood of an arbitrary submanifold embedded in a (pseudo- 
)Riemannian manifold. This can be described by Fermi normal coordinates (FNC) satisfy- 
ing certain conditions as described by Florides and Synge in [7] . By generalizing the work 
of Muller et al in [22] on Riemann normal coordinate expansion, we derive all order FNC 
expansion of vielbein in this neighborhood with closed form expressions for the curvature 

t— i expansion coefficients. Our result is shown to be consistent with certain integral theorem 

for the metric proved in [7]. We explain how this result is relevant in understanding an 
interacting multi-particle system forming bound states in a curved background and a loop 

,-h space description of non-linear sigma model. 
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1 Motivation and summary 

As it is well known, Riemann normal coordinate (RNC) system is a very useful tool in 
differential geometry and general relativity [HE]. In particular, it is used for computing 
covariant Taylor expansion of tensors around a point in a manifold [31 H] • 

Fermi normal coordinate (FNC) system is a generalization of RNC in the sense that 
the RNC-origin, which may be viewed as a zero-dimensional submanifold of the ambient 
space, is replaced by a higher dimensional onq]] It is related to a general coordinate sys- 
tem through the exponential map along the normal directions. The region surrounding 
the submanifold where FNC is well defined, i.e. where the exponential map is a diffeomor- 
phism, is called a tubular neighborhood. The existence theorem for such a neighborhood 
suffices as a powerful tool for various analysis in differential geometry [Sj. 

The use of FNC in general relativity has been mainly due to the fact that one wishes to 
study various physical effects around a particle worldline (one dimensional submanifold) in 
a given background. A partial list of references is [HI QUI IH [111 HE] . In all such applications 
one needs to know explicitly the covariant Taylor expansion coefficients of geometric 



X FNC was first introduced by E. Fermi for a curve in 5j. It was then followed by various other 
generalizations (see [6] and references therein). Our consideration, which is same as that of |7], will be 
described in detail in 331 



quantities around the worldline/submanifolcQ So far only a few such coefficients have 
been know. 

In one line of research [13] FNC expansion around a worldline has been used to com- 
pute curvature corrections to the energy spectrum of hydrogen-like atoms in curved back- 
grounds. The reason for considering a worldline is the following approximation (see also 
[To]): the center of mass (CM) of the system falls freely along a geodesic during the time 
interval relevant for the study of the system (e.g. an atomic transition). This is of course 
a fair assumption when the CM mass is large. However, a framework for a more com- 
plete treatment needs to be investigated. For example, one may wonder how the internal 
dynamics of the atom (or any other composite object) backreacts to the motion of the 
CM. In other words, if it should be possible to derive the geodesic equation for CM at 
the leading order in a semi-classical expansion^ where the CM mass is large. 

To answer the above question one needs a fully relativistic treatment of a multi-particle 
interacting system forming a bound state. This is a difficult problem [Tj5]. However, it is 
interesting to note that in spirit it is very similar to the computation of a' corrections in 
string theory, though there are differences in details (of course). Here the analogue of the 
composite system is a closed string to be viewed as a bound configuration of an infinite 
number of particles. The relativistic theory governing the system is given by the two- 
dimensional non-linear sigma model (NLSM). The analogy can be seen very explicitly 
in a loop space description studied recently in [20j EI]- In particular, in pTJ a semi- 
classical limit {a' — > 0) with the following property has been discussed. If the target 
space is M then the semi-classical expansion is related to the tubular expansion of NLSM 
around the embedding M M- LM, where LM is the ( infinite dimensional) loop space 
and the submanifold, which is isomorphic to M, is the space of all vanishing loops. The 
latter turns out to be the space of all possible positions of CM of the loop. The internal 
dynamics is described in the tubular neighborhood and the effective theory of the CM in 
M is obtained by using a Born-Oppenheimer type approximation. 

Since a' corrections to the effective theory needs to be derived from the tubular expan- 
sion of NLSM, it is important to know the Riemannian structure in the relevant tubular 
neighborhood. It was shown in [21] that such a Riemannian structure is related to Rtm, 



2 The same is true for certain applications outside physics also. See [H] which describes how this is 
used in computing volume of tubes \T7\ [18] . 

3 What we have in mind is an expansion in inverse of the CM-mass. The quantum mechanical h may 
be set to 1. 



i.e. a particular Riemannian structure on TM that views the zero section TMq(= M) as 
an embedded submanifold admitting a tubular neighborhood. 

In £J2]we will discuss how understanding Htm is also a problem relevant to the original 
multi-particle system in M. Moreover, we will describe two steps following which Htm 
may be derived completely. The first step is to have an understanding of tubular expansion 
in a generic case where one considers an arbitrary submanifold embedded in a (pseudo- 
)Riemannian ambient space. The second step is a mathematical problem [23J that we will 
make precise in ^J2] itself. In this work we would like to solve the first step. In particular, 
we will derive all order covariant Taylor expansion of the vielbein components in such a 
generic case. Using this result one can compute the metric-expansion to any desired order. 
In fact, we will present explicit results up to 10-th order. This enables one, in particular, 
to extend the analysis done for the general relativistic applications by considering particle 
worldline as mentioned earlier to higher orders. 

We now briefly discuss the technical points relevant to our analysis towards pursuing 
the first step mentioned above. In [7J, Florides and Synge (FS) constructed the special 
coordinate system under consideration for an arbitrary submanifold embedded in a higher 
dimensional (pseudo-) Riemannian spacaj They wrote down the special coordinate condi- 
tions in terms of the metric and proved an integral theorem describing its behavior away 
from the submanifold. We will review the basic results in £J3j 

In [22], Muller, Schubert and van de Ven considered the RNC coordinate conditions 
written in terms of the vielbein and spin connection. Using certain differential geometric 
techniques it was possible to write down an integral equation for the vielbein in terms of 
the Riemann curvature tensoi|j The authors were able to solve this equation to produce 
the complete RNC expansion of vielbein with a closed form covariant expression for the 
curvature expansion coefficients. 

Here we will use the same techniques to generalize the results of [22] to the case of FNC. 
In particular, we derive the integral equation for the vielbein and closed form covariant 
expressions for its curvature expansion coefficients in the set up considered in [7J. The 
results are different when the vector index of the vielbein takes values along the directions 

4 The authors of [7] called it a submanifold based coordinate system. However, following the modern 
nomenclature (see, for example [TB]) we will continue to call it FNC. 

5 As mentioned in [22], such coordinate conditions are the gravity analogue of the Fock-Schwinger 
gauge [24 in gauge theory which can also be used to write down an integral equation for the gauge 
potential in terms of the field strength [25] . 



tangential and transverse to the submanifold. The transverse results are exactly the same 
as that of [22], as expected. This will be discussed in £J4J In §[5j we show how our result 
is consistent with the metric integral theorem of (7J. Explicit numerical results for the 
expansion of vielbein and metric up to 10-th order in FNC have been presented in one of 
the appendices. 

2 On the relevance of Rtm 

Here we will explain how Htm is relevant to the study of a multi-particle interacting 
system forming a bound state in a curved space. This will also specify Rtm more precisely. 
Then we will explain how the problem of determining Rtm is related to the result of 



tubular expansion (see {4.2) in a generic case found in this work. 

We begin by considering a system of two interacting particles in a curved space, 
say M. Extension of the argument to higher number of particles is straightforward. We 
assume that there exists a generally covariant quantum mechanical theory which efficiently 
describes the bound state and that the configuration spacqj of this theory is a normal 
neighborhood of M x M [] which corresponds to having both the particles in a single 
normal neighborhood of M. The latter implies, in particular, that the diagonal A of 
M x M intersects with the configuration space. 

We now argue that the CM lies in the diagonal submanifold A(= M) 4Mx M. This 
is evidently true when the two particles coincide in M. When the particles are separated, 
one may proceed as follows. Given a two-particle configuration, erect a tangent space 
T X M such that x lies in the neighborhood of that configuration. Then find the preimage 
of the configuration in T X M under the exponential map exp^, : T X M — > M. If the sum of 
the position vectors (in the preimage) in T X M vanishes, then x is the CM. Therefore all 
possible configurations can be described on the tangent bundle TM in the following way. 
The CM always lies on the zero section TM (= M) whereas the relative position is given 
by an element in the fibre positioned at the CM in TMq. Notice that in this description 
any given two-particle configuration lies entirely on the same fibre, which is not true in 



6 We mean extended configuration space, as M includes time. 

7 In the following discussion we will repeatedly encounter such spaces. For simplicity and to avoid 
clutter of notations, we will refrain from emphasizing each time that our considerations are local, i.e. 
restricted to a suitable normal/tubular neighborhood. 



the original description in M x M. The above procedure defines a bundle map, 

<$>:TM -> M xM , (2.1) 

which relates the two descriptions of the configuration space. It is a difleomorphism (given 
that exponential map is a diffeomorphism) and, 

$ : TM -> A . (2.2) 

This shows that A is the space of all possible CM's. 

Given the above discussion, it is now clear that the internal dynamics of the system 
lives in the tubular neighborhood of A c — >■ M x M. In order to derive the effective theory 
of the CM in A(= M) one may follow the procedure adopted in |2T] for NLSM, namely 
one first expands the generally covariant quantum theory in the relevant FNC and then 
uses a Born-Oppenheimer type approximation described in that work. 

We now specify Htm precisely. The FNC relevant for the embedding A c -^ M x M 
is related to the general coordinates on M x M (which gives the usual direct product 



Riemannian structure) by a coordinate transformation which is precisely given by (2.1). 
In other words, the Riemannian structure given by the FNC expansion of the vielbein 
can be thought as the Riemannian structure Htm that has been induced on TM through 
the above procedure. It is such that TMq(= M) is viewed as a submanifold of TM which 
admits a tubular neighborhood/expansion. 

We now specify the mathematical problem of finding Htm completely, given the results 
of this work. To this end we note the following two points, 

• The relevant tubular expansion in TM is a special case of the general form derived 
in this work. As displayed in §4.2[ the expansion coefficients are tensors of the am- 
bient space (in this case TM) evaluated on the submanifold with specific numerical 
prefactors. Such tensors are various combinations of curvature and its covariant 
derivatives. These coefficients carry detailed information about the extrinsic prop- 
erties of the embedding. 

• In the generic case considered here, the intrinsic and extrinsic properties of the em- 
bedding can be chosen somewhat independently. However, the Riemannian struc- 
ture of TM is induced from that of M. In other words, there is only one set of 



independent geometric data involved in the problem, which is that of M. There- 
fore, for the embedding TMq(= M) ■=->■ TM, the tubular expansion coefficients must 
all be expressed in terms of the intrinsic geometric data of M. 



Such relations have been found in [2T] up to quadratic order with an undetermined real pa- 
rameter. Determining Htm completely is equivalent to finding all such relations. Clearly, 
the all-order-result found in this work is necessary for such an exercise. This is a work in 
progress [23] . 

The rest of the paper is devoted to finding all-order-result for the expansion of vielbein 
components around an arbitrary submanifold embedding. 

3 Metric-integral-theorem due to Florides and Synge 

We consider a .D-dimensional submanifold M embedded in a higher dimensional pseudo- 
Riemannian space L of dimension d. Our notation for indices is as follows: Greek indices 
(a, (3, ■ ■ •) run over D dimensions, capital Latin indices (A, B, ■ ■ •) run over (d — D) trans- 
verse dimensions and small Latin indices (a, b, ■ ■ ■), over all (d) dimensions. Other nota- 
tions and conventions essential for the rest of our discussion are presented in appendix 

m 

The work of [7] proved the existence of a coordinate system z a = (x a , y A ) (that will be 
called FNC following [16]), where x a is a general coordinate system in M whose embedding 
is given by, 

y A = , (3.3) 

and the following conditions are satisfied for the transverse coordinates, 

g a B(x,y)y B = r) aB y B , (3.4) 

where g a b(x,y) is the metric tensor in FNC. 

The following integral theorem was proved in [7], 

g A B(x,y) = r] AB + 2y c y D dt F^lACDB^.ty) , 

Jo 

9Ap{x,y) = y c gAp,c(x) + 2y c y D dt F 2 (t)l A cDp(x,ty) , 



-1 



9ais(x,y) = G a/3 (x)+y g a >c (x)+2y y / dt F 3 (t)l aC Dp(x,ty) . (3.5) 

Jo 

For any function f(x,y) in L, we have defined: f(x) = f(x,0). A comma in the suffix 
indicates ordinary derivative with respect to the argument. For example, 

d 
9ap,c{x) = km ^-^gapix, y) . (3.6) 

G al 3 = g a /3 is the induced metric on M and the functions Fi(t) are defined as follows, 

F 1 (t)=t(l-t) , F 2 {t)= l -{l-t 2 ), F 3 (t) = l-t,. (3.7) 



Finally, 



l aC db(x, y) = -(g a b,cd + g C d,ab - g a d,cb - g C b,ad)( x i v) > ( 3 - 8 ) 



is the linear part of the covariant Riemann curvature tensor. 

4 FNC expansion in arbitrary tubular neighbourhood 

The above metric-integral-theorem involves the linear part of the curvature tensor. This 
form is not very useful for deriving the curvature expansion. Below we will use the 
technique of [22] to derive the curvature expansion for vielbein. In §[5] we will show how 
our result is consistent with that of FS. 

4.1 Integral equations for vielbein 



We will first derive the integral equations satisfied by the vielbein (eqs.(4.24) below). 
To this end, we define, following [22J, the radial vector field (see appendix [A| for our 
notations), 



r = y A £ A ■ (4.9) 



Then the coordinate condition (3.4) can be rewritten in terms of the vielbein and the 



connection one-form in the following way, 



ir 0(o) = s (a) B y B , (4.10) 

W a \ b) = 0, (4.11) 



where ir denotes the interior product [26J with respect to the vector field r| 

The next step is to relate the second order Lie derivative of 9^ a ' and the curvature 



two-form. Using eqs. (A. 58, 4.11) one writes, 

C Y 9 {a) = u {a) {b) ir9 (b) + di r ^ a) . 



Then using the following result, which can be easily proved by using (4.10), 

£ r diJ {a) = dij {a) , 



one arrives at, 



£ r (£r - l)0 (a) = C r co ia) (b) i r 9 (b) . 



(4.14) 



(4.15) 



(4.16) 



Finally, one calculates two sides of the above equation independently. The left hand side 
is directly calculated by noting, 



0(°> = e^ a (x, y)dx a + e^ A (x, y)dy / 



(4.17) 



and the right hand side can be calculated by using eqs. (A. 59, 4.11). This leads to the 
following second order differential equation for the vielbein, 



d(d + e)e™ b (x, y) = p w (c) (x, y; y)e^ b (x, y) 



(o), 



,( c ), 



where 



P (a \b){x,y,y) = r ia) C D(b){x,y)y c y D . 
For any function f(x,y) we have defined, 

df(z,y) =y A d A f{x,y) , 



(4.18) 



(4.19) 



(4.20) 



To see how (4.10 4.11) and ( |3.4[ ) are equivalent one first shows that (3.4) implies, 

j a BC (x,y)y B y c = o, 



(4.12) 
and vice versa. One then rewrites the above equations in terms of the vielbein and spin connection 



coefficients to derive (4.11) and, 



a,.b 



d A eW B (x,y)y A y" = 



(4.13) 



which, in turn, is solved by (4.10) 



and, 



1 when b = B , 
■1 when b = (3 . 



(4.21] 



Following [22], we Taylor expand both sides of (4.18) around y = and equate the 
n-th order terms. This gives, 



d£e (a) B (x,0) 



n — 1 



71+1 

d£e%(x,0) = dr 2 [p (a \c)(x,0;y)e {c) p(x,0)] , n>2 



n > 1 



(4.22) 



where we have introduced a new notation, 



d n [f(x,0;y)g(x, 



y 



,A n 



y An lim 



d 



d 



[f(x,y;y)g{x,y)] . (4.23) 



y->6 dy Al dy A 
The integral forms that solve the above equations are as follows, 

e^ B (x,y) = 6&B+ t dtF 1 (t)p^\ c) (x,ty;y)e^ B (x,ty), 

Jo 

e^ {x,y) = e^\{x)+y A e^\ A ^)+ [ dt F 3 (t)p^\ c) (x,ty;y)e^ p (x,ty) . 

Jo 

(4.24) 

The first term in the first equation and the first two terms in the second equation are 
"initial conditions". The former is obtained by contracting the first equation with y B , 



using (4.10) and then noticing that the integral term does not contribute under such a 



contraction because of the anti-symmetry property of the Riemann tensor. For the latter 
we havq^J 



e^f, = E^ p {x) , e {A \ = , y A &\, A = y A 0Jp {a) A , 



(4.25) 



where E^ a K is the vielbein of the induced metric G n p, on the submanifold. While the first 



a/3 



equation is obvious, the second equation is required by consistency with (3.4). The third 



equation is obtained as follows. Using the following relation (total covariant derivative) , 

dj b \ - Ya C e {b) d + tu a {b \ e) e^ c = , (4.26) 



9 In the rest of the discussion in this subsection, we will refrain from explicitly writing down the 
arguments. It is understood that a geometric quantity denoted by a lower case symbol without a bar has 
an argument (x, y) and the same with a bar has an argument (x). 
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and (4.10) one gets on the submanifold, 



~d ^(a) _ -d 2j(a) _ - (a) 



(4.27) 



Then using (4.26), (4.27) and (4.11) one arrives at the third equation in (4.25). 



4.2 Closed form expressions for curvature expansion-coefficients 

We write, 



e^\(x,y) = £eM 6 (x,y) . 

p=0 

where e p is the contribution at p-th order in curvature, such that, 



e o s 



>) 

e /8 



#( a >/j + cV a) c2/ C , for a = a 



W/3 (A) Cl/ C , 



for a = A . 
Following the same techniques as in [22] we find, 

ejf> B (x,y) = J2 ^\si,s 2 ,- ■ ■ ,s p )[(y.V)^p(x,0;y) ■ ■ ■ (y.V) s "p(x,0; 



Sl,—,Sp 



e p a) ^x, y ) 



where, 



J2 jf\si, s 2 , • • • , s p ) [(y.V)Xa;, 0; y) • • • (y.V) s »p(x, 



Sll'",Sp 



(4.28) 



(4.29) 



1 (a) (6) 



(a) (6) 

; y;j W e /3 > 



(4.30) 



J?)( 



q \ s li s 2,---,Sp) 



Si 



-i / 1 rft 1 ^ 1+ - +Sp+2p " 2 F (? (t 1 ) f 1 dt 2 t s 2 2+ - +Sp+2p - 4 F q (t 

■ ' ' s p- Jo Jo 



Q'tptp " g\p; 



VI i 



9=1,3 



g^2j 



(4.31) 



and -F\(£) and F 3 (t) are defined in eqs.(3.7). Explicit calculations give the following results, 

C ( i p \s 1 ,s 2 ,--- ,s p ) 



:(P) 



J 7 ! '(Si, S 2 , ■ ■ ■ , S p ) 



J r i P \si,s 2 ,---,Sp 



[s 1 + s 2 + --- + s p + 2p+l)\ 

C 3 [Si,s 2 , ■ ■ ■ , Sp) 
{s 1 + s 2 + --- + s p + 2p)\ ' 



(4.32) 
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where, 

ri(p) ( „ „ „ N _ ris 1 +s2-\ \-s p +2p-l/^s2+S3-\ \-s p +2p-3 s~is p +l 

^1 \S 1 ,S 2 , • ,S p ) — U Sl " L/ S2 * '"^sp 5 

Ci P \ Sl ,S 2 ,---,S p ) = Q+-2+-+-P+2P-2C--2+-3+-+-P+2P-4 . . . j_ ^ (433) 

C™ being the binomial coefficients. 

In order to present the results in convenient matrix forms we introduce d x d matrices 
Eq, Ej} and R 2+S (x, y) such that their elements are given as follows, 

[AU = eo a) /3 , pEj]aB = , (4.34) 

[E^] a/3 = , [E^] aB = e { a) B , (4.35) 

[R 2+s (x,y)U = (yVyp {a \ b) (x,0;y) . (4.36) 

Notice that R 2 + s is linear in curvature, but (s + 2)-th order in y. Using these we further 
define the following two matrices: 

oo 

E q (x,y) = I + J2 J2 ^ p) (^---,s p )R 2+Sl (x,y)---R 2+Sp (x,y) , g = 1,3, (4.37) 

p=l si,--,s p >0 

where E is the d x d identity matrix. The vielbein matrix is given by, 

E(x,y) = E^x,y)+E ± (x,y) , 
E»(x,y) = E 3 (x,2/)eS, E ± (x, 2 /)=E 1 ( a ;, 2 /)E^, (4.38) 

such that the nonzero elements of E^(x, y) and E ± (x, y) are given by e^ a \(x,y) and 
e^s(x,y) respectively. 
The metric is given by, 

9ab{x,y) = e {a ' ] a {x,y)ri ia i b/) e {b '\{x,y) . (4.39) 

In matrix form we write, 

<?(x, j/) = E T (x, y)r]E(x, y) = # !l + h + /i T + g L , (4.40) 

where 

<?" = (E") T ^Eii = (e5) t g"eS , 

/i = (E'^E^ = (E|J) T iJE^ , 
12 



g 1 = (E ± fr l E ± = «) t GX , 



(4.41] 



are all d x d matrices whose non-zero elements are, 

W l ]ap = g a p{.x } y) , [h] aB = g aB (x,y) , [g^} A B = 9AB(x,y) ■ (4.42) 



Equations (4.38) and (4.41) imply 



where, 



G 11 = EjE 3 , H = EjEx , G 1 - = E^Ex 



E q = V E q . 



(4.43) 



(4.44) 



We will display explicit results for the matrices Ei, E3, G\ H and G 1 - up to 10-th 
order in y in appendix [Bj 

5 Alternative proof of Florides-Synge theorem 



Here we prove the integral theorem in (3.5) starting from our result in (4.18). Analogous 



integral equations for vielbein, namely (4.24), were derived from a set of equations (4.22) 



describing all the necessary transverse derivatives of the vielbein evaluated on the sub- 



manifold. We will first derive the analogue of eqs.(4.22) for the metric (eqs.(5.53) below) 



The integral equations in (3.5) will then follow as solutions to eqs.(5.53) 
Writing g ab = e a .e b = r] {cd) e {c) a e {d \ and using (418) we get, 



d g ab + e de a .e b + e e a .de b = 2(r aCD by y + de a .de fe ) . 



(5.45) 



Next we show that the right hand side receives contribution only from the linear part of 
the curvature tensor, i.e., 



.CD 



C..D 



TaCDbV y + de a .de b = l a cDby V 



(5.46) 



where l abc d is as defined in eq.(3.8). To establish the above equation we first write, 



'abed '•abed "T Qabcd > 

where the quadratic part of the curvature tensor is given by, 

Qabcd = 9ef(l e ad l f bc ~ 7«7m) • 



(5.47) 



(5.46 
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The condition in (4.12) implies that the first term does not contribute in the contraction 
q a CDby C y D - The contribution from the second term can be calculated by using, 



d a 9bcy = (Vbc - 9bc)S a 



c 



(5.49) 



which follows from eq.(3.4). The result is, 

1 cd 
1 

-de a .de 6 



q a cDby c y D 



--g cd (dg ca + g aC 8\ - g cC S° a )(dg db + 



d - 9c1d5 b ) , 



(5.50) 



establishing (5.46). To arrive at the second line we used the following result, 



de a .e b - e a .de b = y c {g bd l d Ca ~ 9adl d 



Cb) i 



c 



9aC$ b ~ 9bC$ a 



(5.51) 



where the first line is obtained by using an analogue of eq.(4.26) with indices properly 



placed and (4.11). The second line results from a direct calculation using the explicit 



form of the Christoflel symbols in terms of metric. 



Using (5.46) and (5.51) in (5.45) we now rewrite the relevant second order differential 



equations for various metric components in the following way, 



d 2 g A B + dg AB 
d 2 g a/3 - dg afS 



2y C y D UcDB , 
2y c y D l ACD p , 

2y c y D LcD? ■ 



(5.52) 



Following the method of £4.1 we Taylor expand each of the above equations around y = 



to find the analogues of eqs.(4.22). The results are as follows, 

1 



do9Ap(x,0) 



.n 

+ 
2n 



d^ gAB (x,0) = 2— y^X ~ 2 Iacdb(x } 0) , n>l 



y c y n d^l ACD p(x,0) , n>2 



n + V 
d%g af} (x,0) = 2y c y D d n -% CD p(x,0) , n>2. 



(5.53) 



The integral equations that solve the above equations are precisely the ones in eqs.(3.5). 
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A Notations and conventions 

We denote the coordinate and non-coordinate bases of the tangent space by £ a and £( a ) 
respectively, 

£ b = e^ b S ia) , (A.54) 

e( a \ being the vielbein-components. We denote by r\ the diagonal matrix with indicators 
of the non-coordinate basis as diagonal elements, 

Vab = (£(a),£(6))- (A.55) 

The bases dual to £ a and £( a ) are denoted by dz a and 9^ respectively We consider the 
torsion- less situation and denote the Christoffel symbols by •y a bc , 



V b S c = V Sb £ c = ^Sa , (A.56) 

is the Levi-Civita cc 
given by, 



where V is the Levi-Civita connection. The connection one-form ui ,-, = W/j i c \0^ a ' is 



V (a) £ {c) = V i(a) S {c) = uj {a f\ c) S {b) . (A.57) 

Cartan's structure equations in our case are given by, 

d§^+u) ia \ b) A§ (b) = 0, (A.58) 

d^\ b) +J a \ c) AJ% = r«,, (A.59) 

where r L = \r iLtM^^ A ^ * s the curvature two-form, 

r% c)id) = (8 {a \ (V (C )V M - V( d) V (c ) - V^^])^)) , (A.60) 

[., .] being the Lie bracket 



B Coefficients up to tenth order 

Here we will display the expansion coefficients up to 10-th order in y for the matrices E l5 
E 3 , G", H and G -1 . We first expand these matrices in the following way, 

oo 

E q = E + ^Ef(x,i/), 9 = 1,3, 

p=2 
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oc 



p=2 p=2 p,q=2 

oo oo oo 

p=2 p=2 P,9=2 

oo oo oo 

G ^- = V + Y^ G^ = V + J2^i ] + E S P)T ) + ]T Ef )T E? , (B.61) 

p=2 p=2 P,<?=2 

where E,q(x,y), G^ p \x,y), H^ p \x,y) and G ±( - P \x,y) are contributions to the relevant 
matrices at the p-th order in y. The results up to p = 10 are given below. The following 
notation will be used in writing down results for G", H and G ± . Given a matrix of the 
form: R m R n ■ ■ • = i]R m R n ■ ■ •, we define, 

{R m R n • • •} = rjR m R n ■ ■ ■ + 77 • • • R n R m . (B.62) 



Results for Ei: 



Ef = ^R 2 , 



6 



Ef } = -R 3 , 
1 12 d ' 

Ei = R.4 4~ Rn , 

1 40 120 2 ' 

11 1 

E| 5) = R 5 + R3R2 + R 2 R 3 , 

1 180 180 360 ' 

Ei = Rk H R4.R9 + R9R4 + Ro H R9 , 

1 1008 6 504 4 2 1680 504 3 5040 2 ' 

11 1 1 1 1 
Ei = — — — R7 H R^Ro H R4R3 H R3R4 4~ — — — R^Rr, + RoRt; 

1 6720 7 2016 5 1344 4 3 2240 3 6720 3 2 10080 5 

H R9R3R2 H~ R9R3 , 

10080 20160 2 

Ei = Ra H RrR2 H R.^Rs H — z — R/i H — ~ R4R9 

1 51840 10368 5184 5 3 5760 4 17280 2 

111 I 

12960 12960 3 25920 72576 

H R2R4R2 4" R2Rq H R0R4 + R9 , 

36288 36288 3 120960 2 362880 2 ' 

Ei = Rq -\ R7R2 H RfiR-? 4~ R5R4 + R5R9 + R4R5 

1 453600 64800 25920 21600 64800 2 32400 

H R4R3R2 ~\~ R4R9RS -|- — HqUfi + .TlsRaR? ~\~ Rq 

32400 64800 90720 45360 45360 3 
151200 3 453600 3 2 604800 7 181440 5 120960 ' 
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I 1 2 1 „ 1 2 

H rtorts-rLi -|- R2RSR9 -|- — — KnKs -|- — — RqRqR2 

201600 604800 2 907200 2 907200 2 

H RnR^ , 

1814400 2 ' 

Ei = Rm H — z RsR? + R7R1 H RkRa + RfRo 

1 4435200 ° 475200 8 158400 7 3 105600 6 316800 6 2 

II 11 

H Rk H~ R^RqRo -|- — — R^RoR^t -|- — — R4RR 

118800 5 118800 5 3 237600 5 3 221760 6 

H R^R? -|- RaRo -|- R4R0R4 -|- R4R0 -|- — RqR7 

110880 4 110880 3 369600 1108800 2 739200 3 7 

H z — R3R5R2 H z R3R4RQ 4~ RoRa 4~ z RoRo 4~ R3R2R5 

221760 147840 246400 3 739200 3 2 1108800 

H RsR2RsR2 ~\~ z R^RoRs ~\~ — z R2R8 

1108800 2217600 2 5702400 

H R2RrR2 ~\~ — z R2R5RS ~\~ R2R/1 

1140480 570240 633600 4 

H R2R4R9 ~\~ R2RSR5 ~\~ R2RqR2 

1900800 2 1425600 1425600 3 

H R2R3R2RS ~l~ z RoRfi -I- R9R4R2 

2851200 7983360 2 3991680 2 



r>2r>2 , R^R -I- I 

3991680 2 3 13305600 2 4 39916800 



Results for E 3 : 



E 3 4) = i-R 4 + i-R2, 
3 24 24 2 

Eq = Re; -|- RqRo -|- R?Rq , 

3 120 5 40 3 120 3 ' 

11 1 11 

Eq = — Rfi -|- R4R2 -|- z R2R4 ~\~ Rq ~\~ Z Ro 1 

3 720 6 120 720 180 3 720 2 ' 

111 1 1 

Eq = R7 H R5R2 H R4R 3 H R^R4 H R 3 R 9 

3 5040 504 5 504 3 1008 1008 2 

1 11 

H R2R5 ~\~ R2RSR2 ~\~ RnRs , 

5040 1680 5040 2 ' 

11 1 11 

Eq = Rs + R«R2 H R5R3 H R/i H R4R0 

3 40320 2688 2016 2688 4 2688 2 

H — z — R3R.R H~ RqR2 + - z — RSR2R3 H R2Rr 

6720 2240 3 6720 40320 

1 1 1 1 1 

H Z R2R4R2 "I - R2Rq "I - R0R4 -|- Rn , 

6720 10080 3 40320 2 40320 2 ' 

11 1 1 1 

Eq = Rq H — — R7R2 H RrRs H R5R4 4~ R5R9 

3 362880 17280 10368 10368 10368 2 
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H — — R4R5 H — — — R4R3R2 4" ~~ — R4R2R3 H R^Rfi 

17280 5760 17280 51840 

H R3R4R2 H Rq + R3R2R4 4" R3R9 

8640 12960 3 51840 51840 2 

H R9R7 4" RoR^Ro 4~ R9R4R3 4" — — z~~ R2R3R4 

362880 7 36288 5 36288 3 72576 3 

I 1 1 

-\- — — — — RoR^Ro 4" RoRc; 4" RnR^R? 

72576 3 2 362880 2 120960 2 3 

H R9R3 , 

362880 2 ' 

E, = Rin H R«R2 H R7R3 H RkR4 H RrR 

3 3628800 129600 64800 51840 51840 2 

H Re; 4" R5R3R2 H R5R2R3 H RaRr 

64800 5 21600 64800 129600 

H R^Ro 4" RaRo 4" R4R0R4 

21600 4 32400 3 129600 

H R4R9 4~ R3R7 H R3R5R2 

129600 2 453600 45360 

H R3R4R3 4" — — — R0.R4 + — — — RoRn 

45360 90720 3 90720 3 2 

H RqRoRc; 4" R3R9R3R9 4" R3R9R3 

453600 3 5 151200 3 3 453600 3 2 3 

H RoRs 4~ R^RrRo H RoRc;R^ 

3628800 8 241920 6 181440 5 3 

H R^R/i 4" R9R4R9 4~ R^R^Re; 

241920 4 241920 2 604800 3 5 

II 1 

H RoRoRo H R9R3R2R3 4" RnRp. 

201600 3 604800 3628800 2 
1 11 
H R0R4R2 "I - — z R9R0. 4" R0R4 

604800 2 907200 2 3 3628800 2 



1 
3628800 



+^^^R|. 



Results for G": 



G 
G 

G 

G 

G 

G 



2 • 



« 2 > = R 

! = ^ ' . - 

= — R4 H — R2R2 1 
12 3 ' 

l(5) = — R 5 4- — {R3R2} , 
60 5 60 1 3 J ' 

= Rfi + 1 R4R9 r H R-?Rq + — R9R9 , 

360 6 360 X J 180 3 3 45 2 ' 

l(7) = -*— R7 4- — {R5R2} 4- — {R4R3} 4- ^^{RsR 2 .} 4- — R2R3R2 , 
2520 315 1 J 504 1 3J 2520 1 2J 420 
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G ll(8) = -^R 8 + ^^{R 6 R 2 } + — ^{R 5 R 3 } + — R4R4 + — {R4R2} 
20160 10080 X J 20160 X J 2016 6720 X 2J 

+ 20i60 R3R5 + Joi60 {R3R3R2} i + iooso 11311211 ^ sseo 112114112 ^ 3i5 R2R2 

G " (9) = 1814p^ + 181440 { ^ } + W^^^ 

+ 18M40 {R4R3R2} + 25920 {R4R2R3} + 90720 {R3R4R2} + 6480 RsR ' 

+ 181440 {R3R2} + 18144 R2R5R2 + 181440 {R2R3R ' } 

G ll(10) = — - R 10 + 37 {R 8 R 2 I + —^— {R 7 R 3 1 + -^^{R 6 R 4 } 

1814400 1814400 X J 453600 X J 129600 X J 

+ 1^400 {R6R2} + 1296TO R5R5 + 3024 {R5RsR2} + 907200^ "^ 

+ 604800 {R4R4R2} + 453600 {R4R3} + 129600 R4R2R4 + 1814400 {R4R2} 

11- 2- 13- 59- 
H j RsRkRq r H" R^RdR^ H jRsR^Ror + R 3 RoRsR2 

18144 1 3 5 J 2835 28350 X 2J 151200 

H z R 3 RoR 3 ~\~ RoRfiR9 ~\~ 1 R2R4R0 r ~\~ R2R3R2 

907200 2 60480 86400 X 2J 100800 3 

59 - 2 

H R9RQR2RQ -f - ^— R9R0 . 

151200 3 3 14175 2 



Results for H: 

H(2) = 3^ 2 ' 



4 
hW = A^H-^R 2 R 2 , 

15 15 

1 - 1 - 5 - 

H {5) = ^R5 + ^ R 3R2 + ^ R 2R3 , 

ii = Rr -|- — R4R2 4~ RsR^ 4~ R2R4 4~ R2R0 1 

420 63 140 140 504 2 ' 

1 1 - 1 19 1 - 7 

R7 H R5R2 + R4RS H R^R4 H R^Ro H- R2R5 

2880 480 192 3 2880 320 2 1440 

1 D D D 1 ' D D D 

H (S) 



H^ = 



H ; — R2RSR2 4" R2R2R 3 1 

1440 2880 

1- 1 - 47 1 - 1 - 

= R« 4~ RrRo 4~ R^R-? 4~ R4R4 H R4R0 

22680 8 2835 6 45360 5 3 630 1512 4 2 

H RqRc; H RqRqRo H RqR9R-3 H RoR« 4~ 

11340 3 5 11340 3 3 45360 3 3 1134 2 6 11340 

H R2R0 4~ — R2R2R4 H R2R0 , 

6480 3 7560 2835 2 ' 

ii = Rq H R7R2 4~ — — — RrRs + R5R4 H RkRo 

201600 604800 5760 3200 604800 2 
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H R4R5 -|- R4R3R2 ~\~ R4R2R3 ~\~ JA,S-H,fi ~\~ R3R4R2 

86400 201600 86400 40320 60480 
1- 1- 17- 3- 37 - 

H JA-slxq -|- R3R2R4 "I - R3R9 -|- R2R7 "I - RqRkRo 

1120 3 2688 120960 2 22400 60480 
I - 23 79 229 

-\ R9R4R3 + R9R3R4 -|- R9R3R.} + lvjlv>JA,e; 

840 3 22400 3 201600 3 2 604800 5 

Q1 1 R7 



340 22400 201600 " 604800 

-— R9R9R3R2 ~\~ l\2lv>lv3 i 

37200 3 604800 2 3 ' 

ij ' = Rin H R.sR.2 H R7R3 + — R,f;R,4 H — RrRo 

1995840 ° 1995840 8 199584 7 3 712800 6 9979200 6 2 
1 1 - 23 13 1 
RcRk H R5R3R2 H RkRoRs H R4R« 4~ R4R4R5 



1 - 1 - 23 13 1 

R5R5 H R5R3R9 H RrRoRs H R4R« H 1 

14256 12320 498960 199584 6336 

R4RQ -|- R4R2R4 ~\~ — z R4R0 -|- R3R7 

199584 3 4989600 9979200 2 1995840 

R3R5R2 ~\~ R3R4R3 -\- — — R3R3R4 + — — R3R3R0 

28512 3 5 399168 9979200 9979200 2 

79 - 1 - 101 - 1 

— — R3R2R5 H R3R2R3R2 + R3R9R3 H — — — R2R8 

997920 11088 1995840 2 57024 8 

R2RrR2 H R2R5R3 H — — R9R/1 + R2R4R0 

199584 b 66528 475200 4 1108800 2 

1 - 257 - 29 - 1 

R9R3RC; -|- — — R2R3R9 -|- R9R3R9R3 + — — — R9R9RR 

4455 3 5 997920 3 199584 2 3 2 3 14784 2 2 6 

H R9R2R4R9 H~ — z — R2R9R-? + — z R9R0R4 -|- R9R9 , 

399168 99792 2 2 3 9979200 2 155925 2 ' 



Results for G^: 

G ±(2) = 1%, 

G^ = %, 
6 



° m = 20^ + Is^ 2 ' 

G ±(5) = 1 



^R5 + ^{R 3 R 2 }, 

1 ^ 17 rT , „ 



yu 40 

G = Rk + 1 R4R2 r H R3R3 H Ro , 

504 6 2520 X J 1008 315 2 ' 

G ±(7) = — R 7 + ^-{R 5 R2} + — {R4R3} + — ^{RsR 2 .} + ^^R2R3R2 
3360 15120 X J 3360 X J 30240 X 2J 15120 

G ±(8) = -^R 8 + ^^{R 6 R 2 } + ^-{R 5 R 3 } + — R4R4+ 113 {R4R2} 
25920 18144 x b J 25920 X J 7200 302400 X 2/ 

1 f\^l 1 1 Q O 

H j R^RqRo r H~ R^RoR^ + R9R4R9 H ^R2Ro 1 

181440 X J 12960 18144 14175 2 ' 

G ±{9) = — - — R 9 + 19 {R 7 R 2 } + — {R 6 R 3 } + -^{R 5 R 4 } 
226800 9 453600 X 7 J 7560 X 6 3J 32400 X 5 J 

H 1 R^Ro f -|- R4R3R0 -|- z - R4R9R3 ~l~ — z — J 

12600 X 5 2J 453600 3 14175 3 90720 



— — — R3 R4 R2 

90720 
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17 - 2 23 1 

H R'hRo H z - RSR2R4 + j R^Ro r H z - R2R5R2 

45360 3 14175 453600 X 2J 5670 

29 - 113 - 41 

+ 90720 R2R4R3 + 4 453600 R2R3R4 + 453600 V**** 1 ® ' 
G±(10) = 2217600 Rl ° + 8553600 {RsR2} + 4435200 {R?R * } + 1108800 {R6R4} + 59875200 ^^ 
+ 356400 R5R5 + 2138400 R5R3R2 + 798336~ R5R2R3 + 6652800 ^ R ' R ^ 

+ 4435200 {R4R ' } + 369600 R4R2R4 + 570240 {R4R2} 

H R3R5R9 H R^R4Rs H — z — z iR^RsRor + z R3R2R5 

2395008 443520 17107200 l 2J 798336 

H z R^R2RsR2 "I - R^RoR^ -I- RoRfiR2 ~\~ R2R. c iR.s 

59875200 19958400 2 855360 2395008 

H z R9R4R0 ~\~ R2R3R5 H~ RoRqR2 ~l" z R2R^R2R^ 

59875200 2 2138400 534600 3 59875200 

H z R2R2R4R2 "I - — R0R0 . 

59875200 467775 2 
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